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Abstract. The collapse of spherical neutron stars is studied in General Relativity. The initial state is a stable 
neutron star to which an inward radial kinetic energy has been added through some velocity profile. For two 
different equations of state and two different shapes of velocity profiles, it is found that neutron stars can collapse 
to black holes for high enough inward velocities, provided that their masses are higher than some minimal value, 
depending on the equation of state. For a polytropic equation of state of the form p = Kp'' , with 7 = 2 it is found 
whereas for a more realistic one (described in Pons et al. (2000)), it reads O.36M0. In 



to be 1.16 (^ 



M^ 



some cases of collapse forming a black hole, part of the matter composing the initial neutron star can be ejected 
through a shock, leaving only a fraction of the initial mass to form a black hole. Therefore, black holes of very 
small masses can be fo rmed and, in par ticular, the mass scaling relation, as a function of initial velocity, takes 
the form discovered by Choptuik (1993)| for critical collapses. 
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1. Introduction 



In the early eighties phapiro and Teukolsky (1980)| ad- 
dressed the question of formation of black holes in as- 
trophysical collapses. Among other, they asked the (the- 
oretical) question of the minimum mass of a black hole 
formed by the adiabatic collapse of a stellar core and "In 
particular, can the effective mass-energy potential barrier 
associated with equilibrium configurations be penetrated 
by low-mass cores with substantial inward, radial kinetic 
energy?". In astrophysical scenarios, black holes can also 
form from accretion- induced collapses of neutron star, if 
the neutron star is part of a binary system or during a 
supernova event, when part of the envelop fails to reach 
escape velocity and falls back onto the new born neutron 
star. Therefore, one can also ask the question: how much 
of inward kinetic energy has to be put to a neutron star 
to make it collapse to a black hole? Can a neutron star al- 
ways collapse to a black hole, provided that it gets enough 
kinetic energy? 

If one looks at static neutron star models, these are sta- 
ble against perturbations if their masses are lower than 
some maximal mass and the density is also lower than 
the "critical" density corresponding to this maximal mass 
(they are located on the so-called stable branch). They can 
collapse to form a black hole when the central density is 
higher than the critical one (see e.g. Gourgoulhon (1991)1). 
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The question is then whether stable neutron stars (i.e. 
with lower density and mass than the critical ones) can 
be "pushed" , with a certain amount of inward kinetic 
energy, to form a black hole, and what is the minimal 
mass of the formed black hole and/or progenitor? The 
problem of the minimal mass of a black hole has been 
solved, first from a more mathematical point of view, by 
[Choptuik (1993) , who discovered the "critical coll apse" 
phenomena in General Relativity (see also Sec. 4.1 for a 
discussion). Still, the problem of a stellar core/neutron 
star with inward kinetic energy have never been com- 
pletely studie d; to our knowledge , only a partial study 
was done by Gourgoulhon (1992)[ The aim of this pa- 
per is to numerically follow the collapse of stable neutron 
stars, with inward velocity profile, and determine the ini- 
tial conditions necessary to obtain a black hole as result 
of the collapse; the purpose is also to allow for a realis- 
tic equation of state (EOS), to study dependence on the 
EOS, as well as on the initial velocity profile, giving ki- 
netic energy. Therefore, emphasis shall be put on neutron 
star properties in General Relativity and the question of 
the physical mechanism giving this kinetic energy will not 
be addressed. 

The complete model in General Relativity is described 
in Sec . |2|, including the system of partial derivative equa- 
tions (2d), the two EOS used (2^), and the procedure for 
obtaining initial numerical models ( |2.3| ). Time evolution 
is studied in Sec. ^, where two different possible evolu- 
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tions for neutron stars are given ( p.l\ and |3.2| ) . Numerical 
results are displayed in Sec. m together with a link with 
the "critical collapse" paradigm ( |4.l| ), properties of initial 
data (4.2) and dependence on the parameters (4.3) are 
discussed. Finally, Sec. summarizes the results and gives 
some concluding remarks. 



2. Evolution of spherically symmetric neutron 
stars 

2.1. Field and Matter equations 

The equations for the evolution of the matter and grav- 
itational fields are derived from the Einstein equations 
using the simplifying assumption of spherical symmetry. 
Space-time is foliated into spacelike hypersurfaces St and 
equation are written in terms of t he classical 3-1-1 forma l- 
ism of General Relativity (see e.g. Arnowitt et al. (1962) ). 
The real parameter t is called the coordinate time and it 
can be shown that, making the choice oi polar time slicing 
and radial gauge, in spherical symmetry the metric can be 
expressed in the form (generalization of the Schwarzschild 
metric, see also Gourgoulhon (1991)): 



-N^{r,t)dt^ 



gf.^dx'^dx" 
-A'^{r,t)dr^ 



r^{dd^ 



(1) 



where N{r,t) is called the lapse function. The two func- 
tions N and A will often be replaced by v{r, t) and m{r, t), 
defined aq^ 



In(A^) and A = 1 - 



2m 



-1/2 



(2) 



The 4- velocity of the fluid is denoted v^ , the fluid radial 
coordinate velocity being thus: 

dr v^ 

and the "physical" fluid radial velocity U, as measured 
locally by the hypersurfaces observer is deflned by: 



Ad. 
N dt 



(3) 



Following Romero et al. (1996), the stress-energy ten- 
sor was taken to be that of a perfect fluid; so hydrody- 
namical equations in General Relativity (V^T/;' = 0) have 
been written in the form of a system of conservation laws 



du 

'dt 



1 d 



,N 



r^^fiu) 



A' 



s{u), 



(4) 



where u = {D,^,t} is the vector of evolved quantities, 
f{u) the vector of fluxes and s{u) the vector of sources 
(for details see Romero et al. (1996)| ). The evolved quan- 
tities are defined from hydrodynamical variables, baryon 



^ We use geometrized units in which the speed of light c and 
Newton's gravitational constant G are equal to unity 



and total energy densities in the fluid frame [hb and e) 
by: 



D = ATuB, 
/i - iE+p)U, 
T = E-D, 



(5) 



where p is the fluid pressure given by the EOS (see 2.2), 
r = (1 — U^)^^/^ is the Lorentz factor of the fluid, and 
E = T^{e+p)-p. 

The gravitational held equations in radial gauge, polar 
slicing and spherical symmetry reduce to two equations 
(with no time evolution of the gravitational field): 



dm 
dr 
dv 

dr 



inr'^E, 



A- 



(^^ + Anrip + (E + p)U^ 



(6) 
(7) 



The system therefore consists of three evolution equations 
for the hydrodynamical variables (Eqs. ^ and two con- 
straint equations for the gravitational field (Eqs. (^)-(0)). 
Finally, an equation of state p{nB, e) closes the system. 

2.2. Equation of state 

In this work two different EOS for nuclear matter have 
been used, to describe microscopic properties of neutron 
stars. Apart from the initial condition models, which were 
calculated with an equation of state p(ti_b) for cold, cat- 
alyzed matter, the used equations of state were of the form 
p{nB, e). The first EOS used is the well-known ideal gas 
model: 



P= [l - l)(e - ruBriB) 



(8) 



where tub = 1-66 x 10^^^ kg is the baryon mass and 7 is 
the index. We chose 7 = 2 which may mimic relatively well 
the properties of neutron star matter but always keeps a 
sound speed lower than c {causal EOS). This EOS has 
been used not only because it is very convenient from a 
numerical point of view, but also because it allows for 
the calibration and test of the code. For obtaining initial 
conditions, the following expression has been used 



p{nB) = KnamB ( 

no 



(9) 



with K = Q.\ (as in the study of rotating neutron star 
models by Salgado et al. (1994) ) and 



TiQ = O.lfm 



-3 



(10) 



To get a more realistic EOS for nuclear matter interac- 



tions, we used the EOS described in Pons et al. (2000). It 
is a relativistic field theory model supplemented by non- 
linear scalar self interactions. Nucleons (n,p) interact via 
the exchange of tr— , uj— and p— mesons. The contribution 
from leptons is given by its non-interacting form, since 
their interactions give negligible contributions. Moreover, 
the star being at chemical equilibrium with respect to the 
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weak processes, the neutrino chemical potentials are zero. 
However, as a difference from Pons et al. (2000) , kaon in- 



teractions are not considered in this work. The EOS was 
tabulated and, during the time integration, the interpola- 
tion has been done using bi-cubic splines. This may not be 
as precise as the method presented by 3westy (1996) (us- 
ing hi- quintic interpolation), but thermodynamical con- 
sistency is still preserved. At low densities, this EOS has 
been smoothly joined with a polytrope. As far as initial 
conditions are concerned, the same model has been used, 
but the temperature has been set to zero. 

2.3. Initial Conditions 

To get initial numerical conditions, we first obtained stable 
spherical neutron star models and then, to add a velocity 
profile. A stable neutron star is defined by the fact that 
its central density is lower than the critical one («-™'), 
defined by {Mg being the total gravitational mass of the 



star and n'v 



dM„ 



dn'i 



nBir = 0)y. 



= 0. 



(11) 



The stable model is easily computed, integrating the 
well known Tolman-Oppenhcimer-Volkoff (TOV) system, 
which is the static (all d/dt terms and U are set to zero) 
limit of momentum evolution equation in system (0) , plus 
the equations for gravitational fields Eqs. (p|)-(R). For the 
polytropic equation of state ['j = 2,K = 0.1), one has 



^crit 



3.18 no and the maximal stable mass is A/g 



3.I6M0. For the tabulated EOS of [Pons et al. (2000)| (see 
previous section), one has n^'* = 10.5 tiq and the maximal 
stable mass is Mg = 2.O8M0. The rather high value for 
maximal mass for the 7 = 2 polytropic EOS is linked to 
an also high value chosen for K. Thus mass scales between 
both EOS used in this work were different. However for 
polytropes of the form (|^) with 7 = 1 + 1/n, the 
constant A'"/'^ has units of length in geometrized 
units. Following Cook et al. (1994), 



one can use 



this constant to set the fundamental length scale 
of the system. Defining dimensionless quantities 
as in [Cook et al. (1994)| (t = R-^'/H, p = A'>, 
ub = K^'ub, ...), one can see that in our case all 
masses (even in dynamical evolution) scale like K^-^ . 
Once the density and gravitational fields are com- 
puted, the inward velocity profile is added in the r.h.s. 
of gravitational field equations (0)^(0). The metric coeffi- 
cients are then relaxed in order to take into account the 
contribution from the kinetic energy in the total gravi- 
tational mass. Therefore, initial conditions are consistent 
with the gravitational field equations. The considered pro- 
files were of the form 



U{r) 



A{t) 



V 



N{r) \i?surfacc, 

where V{x) has one of these forms {x — r/i?surfacc) 



Vix) 



K, 



-(x^-Sa;) 



V{x) 



27V. 



amp 



(13) 



10\/5 V^" 3 
(Kimp is the parameter defining the amplitude of the pro- 
file). 



Both are such that the minimal value in [0,1] is —Va 



amp J 



but (|2D verifies V'{1) = and (|l3|) has a null diver- 
gence at the surface. These profiles correspond to what 
is usually observed in collapses of neutron stars (see 
e.g. Gourgoulhon (1991) ) or in the formation of neutron 



stars (see e.g. Romero et al. (1996) ). The initial mini 



mal value of the velocity U{r) in the interior of the star 
(r e [0, i?surfacc]) IS notcd — C/min- It shall also be called 
the initial velocity profile amplitude. 

3. Dynamical scenarios 

Numerical time integration has been done using 
High Resolution Shock-Capturing schemes (HRSC, see 
Banyuls et al. (1997)) for the hydrodynamical system 
(W). The metric constraint equations (0)^(0) were in- 
tegrated using standard finite-differences methods. The 
need of using numerical methods able to handle shocks 
comes from the fact that strong discontinuities can 
form (see Sec. (3.2 )), as it has first been observed by 
[Gourgoulhon (1992) , Unfortunately, this last study has 
been limited by the use of spectral methods^ unable to han- 
dle shocks. 

3.1. Direct collapse 

For a stable star close to the maximal mass, but with 
an inward velocity profile, there may be two final issues. 
If the velocity is relatively small in amplitude, the star 
enters an (theoretically) infinite series of oscillations, for 
no viscous nor radiative damping is present. With higher 
velocity, a collapse to a black hole occurs "normally" , al- 
most l ike for an unstable n eutron star configuration (as 
in e.g. [Gourgoulhon (1991) ). Figure ^ shows energy (E) 
and velocity profiles at different moments of the col- 
lapse of such a neutron star. It corresponds to an initial 
configuration of 2.74 Mq, where the initial velocity pro- 
file is of the type (pj^), with Famp = 0.4c, which gives 
an initial minimal radial velocity of —C/min = —0.6c. This 
collapse has been followed up to i = 0.48 ms, when then 
central value of the lapse became 5.13 x lO^'^, with a ra- 
dial velocity at the outer edge of -0.99 c. Due to the choice 
of polar slicing which avoids the appearance of singular- 
ities, the horizon of a black hole cannot be numerically 
described. However, the final state is a black hole, but 
from the "frozen star" viewpoint, with a collapse of the 
lapse. These criteria have been used in Sec. ^to determine 
whether the result of the collapse was a black hole or not, 
namely a central value of the lapse N{r = 0) < 10~^ with 
an ingoing radial velocity. If both of these conditions are 
fulfilled, it seems unlikely that the collapse lead to any- 
thing else than a black hole^ 



(12) 



Still, a difference from [Gourgoulhon (1991)| in the dy- 
namical evolution can be seen in Fig. [l|: a strong gradient 
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Density of energy seen by eulerian observer 



Density of energy seen by eulerian observer 




rodius [l<m] 
Radial velocity seen by eulerian observer 



radius [km] 
Radial velocity seen by eulerian observer 




radius [km] 

Fig. 1. Snapshots of eulerian energy density E and 

radial velocity profiles U, for a 2.74 M© neutron star de- 
scribed by a 7 = 2 polytropic EOS, with an initial veloc- 
ity profile of type (|l3[), and a starting minimal value of 
the velocity of — C/min = —0.6c. Labels correspond to 
the following times: ti ~ 0.054 ms, ^2 — 0.088 ms, 
is = 0.104 ms, ti = 0.i26 ms, tg = 0.i49 ms, te = 0.202 
ms, tr = 0.225ms, tg = 0.248 ms and tg = 0.271 ms. 
All profiles have been cut at the surface of the 
star, defined by the surface density of the initial 
configuration. 



appears in the velocity profile around r — 9 km, even 
quite before the "frozen star" regime (where the metric 
potentials A and N also exhibit strong gradients). This 
had already been observed by Gourgoulhon (1992)[ where 
it limited the study. It is shown in Sec. 3.2 that, with 



higher initial velocities, a shock can form and part of the 
infalling matter can be ejected. In the present "direct col- 
lapse" regime all the matter ends in the black hole so 
that the mass of the resulting object is that of the initial 
neutron star. In Sec. ^, this has been used to distinguish 
between both possible types of collapses. 




radius [km] 

Fig. 2. Eulerian energy density and radial velocity pro- 
files, for a polytropic •y — 2 1.91 M© neutron star col- 
lapse, with an initial velocity profile of type ([l3|), and 
a starting minimal value of the velocity of —0.62c. The 
displayed profiles correspond to the following times: 
ii = 0.064 ms, ij = 0.123 ms, tg = 0.143 ms, t^ = 0.16 
ms, ^5 = 0.178 ms, tg = 0.216 ms, ty = 0.265 ms, 
tg — 0.348 ms and tg — 0.477 ms. All profiles have 
been cut at the surface of the star, defined by the 
surface density of the initial configuration. 

3.2. Shock and bounce 

For some initial velocities, the "strong gradient" of pre- 
vious section turns into a shock, and part of the infalling 
material is ejected. This can be seen on Fig. where en- 
ergy density and radial velocity profiles are displayed, for 
an initial neutron star of 1.91 Mq (the mass is lower than 
that of Fig. n due to a lower initial central baryon density) . 
The shock appears around t = t^, r — 10 km and then 
moves out of the numerical grid (t^) due to the accumu- 
lation of matter still falling at r > 10 km. The velocity on 
the left side of the discontinuity grows and reaches positive 
values. Later during the collapse, one can see that part of 
the falling matter starts moving outward and reaches es- 
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cape velocity from the central object t^. This of course 
reduces the possibility of forming a black hole, since mat- 
ter is spatially less concentrated. Still, the dynamics of 
the central region (where the velocity is still inward) can 
proceed to a black hole, which shall then accrete matter 
which has not reached escape velocity. This has been ob- 
served in the collapse displayed on Figs. S where, at the 
end of the computation (t = 0.63 ms), the central lapse 
was N{r = 0) = 1.18 x 10^^ and the velocity on the edge 
of the "frozen star" equal to -0.999 c. The central region 
would collapse "directly", as described before in Sec. 3.1. 
The final mass of the black hole was 1.61 M0 (and the 
areal radius in RGPS coordinates equal to 4.77 km), which 
shows that almost 16% of the initial matter of the neutron 
star has been ejected. 

4. Numerical results 

The dynamical fate of neutron stars with initial velocity 
are now studied, keeping in mind the different scenarios 
described in previous section and varying the mass of the 
progenitor, as well as the amplitude of the initial veloc- 
ity profile. This exploration of the parameter space shows 
the possibility of formation of very low mass black holes, 
which can be seen as a feature of the "critical collapse" 
phenomena. 



4.1. Link with critical collapses 

The critical collapse phenomenon (for an interesting re- 
Gundlach (1999)| ) was discovered in the early 90s 
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Fig. 3. Masses of the black holes (+) formed through 
velocity-induced collapses, for a polytropic EOS and with 
an initial velocity profile of type (113) . The initial central 



view see 



by Choptuik (1993)|, who was numerically studying the 



gravitational collapse of spherically symmetric massless 
scalar field. Depending on some parameter of the initial 
conditions (which is generically noted p), the final result 
of the collapse would be a black hole {p large) or the dis- 
persion of the field {p small). He discovered that, when 
fine-tuning this parameter, he could get black holes of ar- 
bitrarily small masses. Moreover, the relation giving the 
mass of the resulting black hole, close to p* the minimal 
value to form a black hole, appeared to be universal in the 
form: 



M 



BH 



c{p~p.r 



(14) 



where C is a constant and a is called the critical exponent. 
The space-time obtained with p = p* shows the very in- 
teresting geometrical property of self-similarity. 

There have been many works on critical collapses since 
this pioneering one and, in particular, the study concern- 
ing perfect fluids by Ncilscn and Choptuik (2000) is of 



particular interest for our study. It is shown that per- 
fect fluid collapses, in ultra-relativistic regime, also exhibit 
"critical" behavior. Our work can easily be connected to 
this one, if one considers, at fixed central density of the 
initial neutron star, the parameter p to be e.g. the param- 
eter T4mp or the amplitude of the initial velocity profile 
Unun- One can then see that, by fine-tuning the parameter 
V'amp (C^nin being a monotonic and continuous function of 



density is O.lno (see Sec. 2.2) for all the runs. These masses 
have been fitted by a formula of type m4) (solid line) . 



Vkmp) of (|l^ ) or ([l3|), one could get black holes of arbitrar- 
ily small masses, as a result of the velocity-induced col- 
lapse. We have not investigated the domain of very small 
black hole masses, for our code was not designed for it. 
Still, starting with a 1.16Mq neutron star and an initial 
velocity profile given by ( [l3|) (and T4mp = 0.79732), the re- 
sult of the evolution was a 3.7 x 1O~^M0 black hole. This 
corresponds to the lowest black hole mass range (about 
1O~^M0) obtained with our code. In order to test the 
code against the relation ([l4|), the central density of the 
initial configuration has been kept constant, and only the 
parameter Vamp of ( p^ has been varied (and thus f/min)- 
The results displayed on Fig. ^ correspond to the 
masses of black holes resulting from collapses, if the cen- 
tral density of the initial configurations is held fixed to 
O.lriQ. These masses are drawn (crosses) as a function of 
the minimal value of the velocity U inside the neutron star. 
The relation (|lj) is verified, at least for "small" masses 
and the best fit (shown on the figure in solid line) cor- 
responds to a = 0.52, C = 1O.4M0 and p, = 0.8705c 
(f^min = p)- These values has shown to depend on the 
particular central density used for the initial configura- 
tions. The fit is poor for large masses and, in particu- 
lar, the relation (nj) breaks down completely near the 
maximal mass for stable static neutron stars (when cen- 
tral density reaches n^'*). The mass scaling expo- 
nent a — 0.52 differs from the result obtained by 

1) who used an 



Neilsen and Choptuik (2000)| {a 
ultrarelati vistic EOS with 7 = 2. Althou gh the re- 
sult from Neilsen and Choptuik (2000) is rather 
universal, the study here starts from very different 
initial conditions (neutron star/exponential distri- 
bution of matter in their case). In particular, here. 
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the maximal mass for stable neutron stars gives a mass 
limit to the problem. Another difference is that the 
parameter used to get the relation (14) is the am- 
plitude of a velocity profile, which has not been 
studied by [Neilsen and Choptuik (200d)|. Finally, 



it has to be stressed out that the family of initial 
conditions used for this study depends on two pa- 
rameters (velocity amplitude and central density), 
the work of ] 



Neilsen and Choptuik (2000) 



assumed 

only a one-parameter family of initial data. The 
role played by the central density as a parame- 
ter is not that of a "critical" one: if one considers 
the family of initial data given by static neutron 
stars (without any velocity profile) of increasing 
central densities, then the value ng'* defined by 
(11) would seem to be the same as p* in control- 
ling the formation of a black hole, but there is no 
critical behavior at this point. Still it is interesting to 
note that the mass scaling relation is valid for velocity 
induced neutron star collapses, at least in the limit of low 
masses (and therefore relatively low central densities). 



4.2. Masses of progenitors - Initial velocities 

Once the EOS and type of velocity profile are chosen, a 
neutron star model serving as initial condition for velocity- 
induced collapse is completely specified by two parame- 
ters: the central baryonic density rig and the parameter 
Vamp or, equivalently, the neutron star gravitational mass 
and the amplitude of the initial velocity profile C/min- For 
each point {rig, C/min) in this parameter space the ques- 
tion of knowing whether the corresponding neutron star 
would collapse to a black hole have been addressed. First, 
it has been found that there exists a maximal value for 
Umin, for which an initial model of neutron star could be 
computed. For higher velocities, it appeared to be 
numerically impossible to compute the metric co- 
efficients taking into account this large amount of 
kinetic energy: the relaxation would fail to con- 
verge, leading to velocity U > c. This may come 
from the way these initial conditions are set: the 
addition of a velocity profile to a static model does 
not give a well defined result. Matter is not at equi- 
librium, whereas the metric potentials are static. 
Therefore an improvement of this study would be 
to consider more realistic initial condition s (as for 
example in Hawley and Choptuik (2000) ), where 
the velocity profile is not set ad hoc but comes from 
a dynamical interaction. 

This maximal value of C/min depends on the central 
density of the neutron star, as displayed on Fig. |4[ one 
can see that for central densities going to zero (and there- 
fore also masses), this maximal value goes to one. This 
figure also shows the various fates of a neutron star in 
the (wr, C/min) plane, for the (7 = 2) polytropic EOS of 
Sec. 2^ and an initial velocity profile of type (|l3|). Curves 
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Fig. 4. Parts of the parameter space plane {rig, C/min) 
defining the fate of a neutron star with a polytropic EOS 
(7 = 2) and an initial velocity profile defined by (|l^). 

ure. A neutron star set to a point of this parameter space 
can become a black hole or not. Still, some region is not 
really interesting to study. The second curve to be dis- 
played selects neutron stars with a mass lower than the 
maximal one, as defined in Sec. 2^. Studying collapses for 
neutron stars with higher masses is not relevant since one 
knows that in that case the neutron star shall end as a 
black hole; the additional velocity profile of the neutron 
star only gives kinetic energy so that the total mass of 
the star become larger than the maximal one for stars at 
rest. It was then observed that this collapse generally pro- 
ceeds "directly", i.e. no shock is present. One more curve 
gives the limit between the direct collapse domain and the 
part of the parameter space for which the star undergoes 
a shock and bounce, as described in Sec. 3.2. Finally, if 



the parameters of the initial neutron star model are below 
the last curve, the black hole does not form, the matter 
being ejected to infinity. The same regions are displayed 
in the {Mg, C/min) plane, on Fig. [|. In particular, one can 
see that there exists a minimal mass for the initial neutron 
star, to form a black hole. This mass is, in the case of a 
polytrope 7 = 2 and an initial velocity profile (13): 



Af ,:„ ^ 1.164M 



o- 



(15) 



are described from top-right to bottom-left part of the fig- 



This figure also shows the minimal neutron star mass for 
which direct collapse (i.e. no ejection of matter) can occur, 
it found to be 2.276 Mq. 

4.3. Dependence on EOS and velocity profiles 

Although the velocity profile ( |l3|) is supposed to occur 
in most of physical scenarios of gravitational collapse of 
compact stars, one has to check the dependence of results 
of previous sections on a change of this velocity profile. 
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Fig. 5. Parts of the parameter space plane (n^, Mg) defin- 
ing tlie fate of a neutron star with a polytropic EOS 
(7 = 2) and an initial velocity profile defined by (|l^). 



Therefore the study of Sec. 4.2 has been undertaken with 
the initial velocity profile for neutron stars being given 
by (|2|). For initial configurations, one can see that this 
change in the velocity profile induces, for given values of 
rig and C/min, a change in the distribution of kinetic en- 
ergy in the star and thus a change in the total gravitational 
mass Alg, with respect to the velocity distribution given 
by (|lj) . Still, qualitative results concerning collapses were 
retrieved: direct collapse or bounce as dynamical scenar- 
ios and formation of low-mass black hole by fine-tuning 
the parameter Umin- The scaling relation (JlJ) has been 
studied for a fixed central density (chosen to be O.lno as 
in Sec. 4.1) and the best fit of results gives the relation: 



M 



BH 



6.98 X (C/min - 0.8991)"''^Mo, 



for which the exponent {a ~ 0-52) is the same as that 
obtained using velocity profile ([l3|). The two other co- 
efficients (noted C and p^. in Sec. 4.1) are different for 
the reasons stated above, concerning the difference in ki- 
netic energy distribution. The borders of different regions 
of Fig. 1^ also slightly change, nevertheless we get as mini- 
mal mass for the initial neutron star to form a black hole 
for a 7 = 2 polytrope and an initial velocity profile dl2) : 



Mr, 



1.155Mp 



0- 



(16) 



This is less than 1% different from the minimal mass for 



velocity profile (13), that is within global numerical error. 
One then may suppose that this quantity is little depen- 
dent on the particular type of velocity profile chosen. 

The second global "parameter" one would like to 
change is the EOS for nuclear matter, since a rather im- 
portant uncertainty exists on properties of neutron star 
matter. In particular, the 7 = 2 polytrope used in this 
study had the advantage of rapid computation as well as 



numerical stability, combined with a rather good descrip- 
tion of neutron matter properties. But if one wants to go 
further, it is then necessary to use a more realistic EOS, as 
for example the one described in Sec. 2.2. As stated in that 
section, maximal mass for static neutron stars is found to 
be M"'* = 2.08Mq, for a corresponding central density of 
rig — 10.5no. Therefore, results from computations 
using this second EOS would show lower masses 
and higher central densities than those using the 
analytical EOS. The parameter space {rig, C/min) has 
been studied, using the velocity profile given by (13). As 
for the 7 = 2 polytrope, neutron stars with realistic EOS 
would collapse to form a black hole, provided that the 
amplitude of the velocity profile be large enough. This 
collapse could occur directly, as in the case of unstable 
neutron stars or with the bounce of a part of the mat- 
ter, which would be ejected, allowing only for the central 
region to form a black hole. The mass scaling relation is 
recovered, for relatively small initial central densities; for 
example, with Ug = 1.5no, the following relation is found: 

Mgg ^ 5.79 X (C/„,i„ - O.7519)O-^iM0, 

which has been tested for black holes in the mass range 
5x 10^'^ — O.TMq. The global topology of regions described 
in figures H and is kept and the minimal mass for a 
neutron star (with matter described by this realistic EOS) 
to form a black hole is found to be 



M„,in = 0.36Mp 



0- 



(17) 



During the collapses leading to very low-mass black holes 
(a few IQ-^Mq), the central regions of the star would 
reach very high densities and it has been checked that 
sound velocity had never become higher than c. 

5. Summary and conclusions 

With an important amount of inward kinetic energy, sta- 
ble neutron stars may collapse to a black hole, overcom- 
ing the potential barrier which separates both types of 
objects. However, there seem to exist an absolute lower 
mass limit, depending on the particular equation of state 
used to describe neutron star matter properties, below 
which neutron stars cannot pass this barrier. This lower- 
limit mass has been found to be ~ 1.16Mq for a 7 = 2 
polytropic EOS, and ~ O.36M0 for a more realistic model 
of nuclear matter. It seems also that these values are inde- 
pendent from the particular initial velocity profile added 
to static neutron star models. But, as it has been noted 



in Sec. 2.3, for the polytropic EOS, like all masses this 
lower- limit mass scales like K'^-^. 

For both studied EOS, in the case when the star col- 
lapses to a black hole it may either proceed directly, as if it 
were an unstable neutron star, all the matter ending in the 
black hole; or there may appear a shock and bounce, eject- 
ing a part of the matter to infinity, so that only a fraction 
of the initial neutron star forms a black hole. With such a 
mechanism, the resulting black hole mass can be arbitrar- 
ily small if one is fine-tuning the amount of kinetic energy 



J. Novak: Velocity-induced collapses 



added to the initial neutron star. This result is in accor- 



dance with works by Choptuik (1993) on critical collapses, 
more precisely the mass-scaling relation (H) also applies 
here, at least for black hole masses not ton rinse tn the 



Gourgoulhon E., 1991, A& A 252, 651 

Gourgoulhon E. 1992, Class. Quantum Grav., 9, S117 

Gundlach C. 1999, Living Reviews in Relativity 1999- 



4, pttp: //www. livingreviews.org/ Articles/ Volunie2/ 1999- 
4gundlach/] 



neutron star critical mass The mass sraling PYpnripnt 
a found for the polytropic EOS is rather different 
from that found, with quite different initial con- 
ditions, by Neilsen and Choptuik (2000) for ultra- 



relativistic 7 = 2 EOS. 

Nevertheless, the velocity necessary to achieve the col- 
lapse to a black hole of a typical 1.4M0 neutron star is 
enormous (see Fig. pi). Even in the case of a smaller K, 



leading to a smaller maximal mass (see Sec. 2.3), this 
velocity always remains larger than 0.1c. It is therefore 
difficult to imagine any physical process which could in- 
ject such an amount of kinetic energy to a neutron star. 
Tidal (anisotropic) effects when the neutron star is in a 
binary system with another compact object may result 
in the destruction of the star (compressed to a "pan- 
cake" , as ordinary stars passing near a black hole in 
Marck et al. (1996)) or in its oscillations. I t has been ar- 
gued (e.g. by Mathews and Wilson (2000) ) that in a bi- 
nary neutron star system, tidal effects could induce such 
a compression (i.e. increase of central density) that the 
stars could collapse to black holes before their merging. 
This is very difficult to achieve, since in our calculations 
it happened that during the collapse the central density 
could reach a value substantially higher than the critical 
density ng'* (in some cases several times rfs'^), whereas 
the star would not end in a black hole, but would rather 
be dispersed or enter an infinite series of oscillations. Even 
in the case of the neutron star instability transition in 



tensor-scalar theory (see Novak (1998)), the velocity is 
never larger than a few percents of c, therefore quite far 
from values displayed on Fig . |5|. Finally, let us mention the 
recent interesting work by Hawley and Choptuik (2000) 
who studied the critical collapses of boson stars, where 
the interaction with a massless real scalar field can result 
in a significant transfer of energy from the field to the star, 
allowing for the collapse to a black hole of stable configu- 
rations. In the case of neutron stars the interaction with a 
gravitational wave seems unlikely to transfer enough en- 
ergy to lead to the collapse to a black hole. 
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